Slow Coarsening in an Ising Chain with Competing Interactions 
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We investigate the zero-temperature coarsening dynamics of a chain of Ising spins with a nearest- 
neighbor ferromagnetic and an n th -neighbor antiferromagnetic interactions. For sufficiently large 
antiferromagnetic interaction, the ground state consists of n consecutive up spins followed by n 
down spins, etc. We show that the asymptotic coarsening into this ground state is governed by a 
multispecies reactive gas of elementary excitations. The basic elementary excitations are identified 
and each decays at a different power-law rate in time. The dominant excitations are domains of 
n + 1 spins which diffuse freely and disappear through processes which are effectively governed by 
(n + I)-particle annihilation. This implies that the ground state is approached slowly with time, as 
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I. INTRODUCTION 



Ising models with nearest-neighbor ferromagnetic and 
more distant antiferromagnetic interactions exhibit rich 
magnetic ordering The competition between fer- 

romagnetism and longer-range antiferromagnetism leads 
to different ordered states and an associated sequence 
of phase transitions as a function of these two interac- 
tion strengths. Such models were originally formulated 
to help describe the complex magnetism of the rare earths 
[^J. Their unusual magnetic ordering is believed to arise 
from the RKKY interaction ||, in which the exchange 
interaction between localized magnetic moments oscil- 
lates between ferromagnetic and antiferromagnetic as a 
function of their separation. Some of the essential con- 
sequences of this complex situation seem to be captured 
by Ising models with competing interactions. 

One of the simplest versions is the axial next-nearest- 
neighbor Ising (ANNNI) model in which there is an 
isotropic nearest-neighbor ferromagnetic interaction and 
a next-nearest-neighbor antiferromagnetic interaction 
along a single axis H. Even in one dimension, this sys- 
tem exhibits non-trivial magnetic properties. For a weak 
antiferromagnetic interaction, the ground state is ferro- 
magnetic, while for a strong antiferromagnetic interac- 
tion there is an "antiphase" ground state which consists 
of 2 spins up, followed by 2 spins down, etc. For a spe- 
cific ratio of these two interactions, an infinitely degen- 
erate ground state arises in which each spin domain is of 
length 2 or greater jfj). 

Given the very different natures of these ground states, 
one might expect that dynamical behavior is also strongly 
affected by such competing interactions. Our goal is to 
understand the kinetics of an Ising chain with nearest- 
neighbor ferromagnetic interaction J\ and n th -neighbor 
antiferromagnetic interaction — J n , when the system is 
also endowed with single spin-flip Glauber dynamics. 
The Hamiltonian of the system is 



For J n > J\/n, the ground state is a sequence of n con- 
secutive up spins followed by n down spins, etc. Our basic 
result is that for J n > J\ the asymptotic approach of the 
system to this alternating ground state can be described 
in terms of a reactive gas of elementary excitations. The 
rate limiting step of the coarsening is governed by the 
(ri + l)-particlc annihilation of the dominant excitations, 
which implies that the system coarsens in time as t~ x / n . 

To provide a context for this work, let us recall the 
well-known example of zero-temperature (T = 0) coars- 
ening in the Ising-Glauber chain with nearest-neighbor 
ferromagnetic interactions only J?j]. Spin flips inside fer- 
romagnetic domains are forbidden at T = because they 
cost energy, while spins at the domain interfaces can flip 
freely, as indicated by 



+ - 



h + + + + ++ , 

since no energy cost is involved. This spin flip is equiva- 
lent to the hopping of a domain wall "excitation" which 
lies between the neighboring misaligned spins. Two do- 
main walls can annihilate when they meet, a process 
which leads to a reduction of the energy in the system. 
In terms of the spins, the annihilation is equivalent to a 
domain which shrinks to zero size via the process 



+ + + + + 



+ + + 



n 



J, 



E 



SiS 



z+n ■ 



(i) 



This correspondence between the ferromagnetic Ising- 
Glauber chain and a gas of domain walls which undergo 
nearest-neighbor hopping and single-species two-particle 
annihilation provides a simple way to understand the 
i -1 / 2 coarsening dynamics of the spin system ||. In the 
following sections we show that the Ising-Glauber chain 
with competing interactions can be understood through 
a similar picture of reactive excitations. 
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II. SECOND-NEIGHBOR INTERACTION 

A. Strong Antiferromagnetism: J2 > Ji 

Let us now study the Ising-Glauber chain with com- 
peting interactions at T = 0. To be concrete, consider 
first the case n = 2, i.e., a near-neighbor ferromagnetic 
interaction J\ and a second-neighbor antiferromagnetic 
interaction — J2. We focus on the case of strong antifer- 
romagnetic interaction, J2 > Ji, as this is the case which 
leads to interesting dynamics. Starting from an initial 
ferromagnetically ordered up state for simplicity, this sys- 
tem evolves to the . . . +H H ••• ground state by 

a two-stage process. First, there is an initial nucleation 
of down domains within regions of up spins. Within a 
ferromagnetically ordered region, there is an energy loss 
AE = 4Ji — 4J2 when a single spin flips to create an 
isolated down spin. This same energy loss arises for any 
nucleation event which occurs within a domain of length 
> 5 when the flipped spin is at least two lattice spac- 
ings from both domain walls. After this nucleation, the 
single-spin domain can grow to length 2 with an energy 
loss AE = — 4 J2 as long as the length of the neighboring 
domain into which this inflation occurs is greater than 
2. At the end of this nucleation stage, therefore, the sys- 
tem consists of ordered ... + H h + • • • regions, as 

well as domains of size 1, 3, and 4. Further, domains of 
size unity exist only if both neighboring domains are of 
length 2. 

These remaining domains now undergo a sequence of 
reactions which ultimately leads to the system reaching 
the ground state. To determine this evolution, first con- 
sider an isolated 3-domain within an otherwise stable ar- 
ray to 2-domains. Since the spin on either edge of the 
3-domain can flip without energy cost, the 3-domain can 
hop isotropically by two lattice sites, as indicated by 

+ + -- + + + + -- 

+ + -- + + +-- + + -- . 

Similarly, an isolated 1-domain within a sea of 2-domains 
also diffuses freely, since either neighboring spin of the 1- 
domain can flip with no energy cost, as indicated by 

--++-++--++-- 

h H 1 h H . 

These freely diffusing 1- and 3-domains are the elemen- 
tary excitations of the system. 

To appreciate the consequences of this statement, con- 
sider an isolated 4-domain within an otherwise stable ar- 
ray of 2-domains. Since there is an energy cost associated 
with flipping either spin in the interior of a 4-domain, this 
process does not occur at zero temperature. If, however, 
either spin at the end of the 4-domain flips, for example, 



-- + + + + -- + + 

hH h + H h+ , 

then the configuration becomes two adjacent 3-domains 
within the stable sea of 2-domains. Each of these 3 do- 
mains can then diffuse freely; the first step of this process 
is indicated by 

-- + + + + + -- + + 

-- + + +-- + + + -- + + . 

One can thus regard an isolated 4-domain as a [3, 3] reso- 
nant state which is formed whenever two 3-domains col- 
lide. This resonance is short-lived, however, since its 
binding energy is zero. 

Continuing this reasoning, consider the evolution of a 
4, 3 pair within a stable sea of 2's. Since both the 3 and 
4 diffuse freely, these two excitations could move apart 
with zero energy cost. On the other hand, the 3-domain 
can shrink and the 4-domain can grow to size 5, also with 
zero energy cost. If this occurs, the 5 is unstable to the 
nucleation event: 5 — ► 2 + 1 + 2. Due to the fact that 
the 4 can be viewed as a [3, 3] resonance, the penultimate 
stoichiometry of this process is 3 + 3 + 3 — > 1, i.e., 3's 
annihilate through triple collisions. 

Finally consider how isolated l's evolve. Since l's dif- 
fuse freely they react only upon meeting another 1 or a 
3. In the former collision, a 4 is formed, 

+ + -- + + - + -- + + -- 

+ + -- + + + + -- + + -- , 

with an associated energy loss AE = — 4Ji. Because the 
4 is equivalent to a pair of 3's, the penultimate stoichiom- 
etry of this process isl + 1 — >3 + 3. Finally if a 1 meets 
a 3, they react to form a stable pair of 2's, 

+ + -- + + + -- 

+ + -- + + -- + + -- . 

with an associated energy loss AE = — 4J 2 - This can be 
viewed as the two-species annihilation process 1 + 3^0, 
since the pair of 2's that are formed become part of the 
stable ground state. 

The underlying stoichiometry of these processes can 
therefore be summarized by the three reactions 

3 + 3 + 3^1 AE= +AJi-U 2 , 

1 + 1 -^3 + 3 AE= -4Ji, (2) 

1 + 3 ^ AE= -4J 2 . 

Since each of these processes leads to a loss of energy, 
they each occur at the same rate when T = 0. The over- 
all effect of these reactions is that the density of these 
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elementary excitations ultimately vanish. Let us now de- 
termine the time dependence for the densities of these 
excitations in the mean- field limit. Using A to denote 
both domains of length 3 and their density, and similarly 
using B for l's, the rate equations associated with the 
reaction scheme (Es) are 



A = -3A 3 + 2B 2 - AB, 
B = A 3 - 2B 2 - AB. 



(3) 
(4) 



A naive qualitative analysis of these equations indicates 
that A's are asymptotically dominant. Thus using A ^> 
B, the above rate equations simplify to 



A ~ -3A 3 - AB, 
B ~ A 3 - AB. 



(5) 
(6) 



Subtracting Eq. (g) from (g) gives A — B ~ — 4A 3 . 
Given A ^> B, we neglect B to give the closed equation 
A ~ -4A 3 . Comparing with Eq. (§) yields A 3 ~ AB. 
This implies B ~ A 2 , which agrees with B ^ A. The 
end results of these considerations are 



A(t )ei -±= and fl^.i, 



(7) 



which are confirmed by numerical integration of the rate 
equations. 

Let us now adapt the rate equations to the case of one 
dimension. Quite generally, we write 



.4: 

B 



-3R 
R 



(8) 
(9) 



Here R is the rate at which A's disappear due to the triple 
collisions and r is the rate at which both A's and B's dis- 
appear due to their mutual annihilation. These equations 
are at the same level of approximation as Eqs. @ and (|^), 
as we again neglect the effect of interactions between two 
B's. The rate R can be shown to scale as A 3 / ln(l/A) |j] . 
The cubic term is just the mean-field rate of triple colli- 
sions between A's, and the logarithmic correction arises 
because one dimension is the critical dimension for the 
three-particle reaction-diffusion processes JlO| . This re- 
action rate gives the rate equation A ~ — A 3 /ln(l/A) 
for three-particle annihilation. This gives the asymp- 
totic behavior for the density A ~ \J\ntjt, which agrees 
with previous numerical and theoretical treatments [p"l| . 
The two-species annihilation rate r can be estimated as 
r ~ B/t, where r is the reaction time for an A and a B 
to meet by diffusion. This reaction time is proportional 
to the square of the distance between a B and its nearest 
A; therefore, r ~ 1/A 2 which gives r ~ BA 2 . 

Subtracting Eq. @ from (g) gives A—B ~ -AR. Since 
we again anticipate that A 3> B, we ignore B to obtain 
A ~ — AR. This relation implies that triple collisions are 
the dominant kinetic mechanism for elimination of A's so 
that this density decays as in three-particle annihilation. 



The factor of 4 in A ~ — AR indicates that four A's even- 
tually disappear after a triple collision - three particles 
are eliminated in the process A + A + A — > B, and then 
the newly-formed B will eliminate another A. To deter- 
mine B(t), note that Eq. (@), together with A ~ — AR, 
imply that r ~ R; that is, the gain and loss terms in 
Eq. (JsJ) cancel. The relation r ~ R can be rewritten as 
B ~ R/A 2 ~ -A/A 2 . We therefore conclude that the 
density of elementary excitations are 
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(10) 
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FIG. 1. Density of domains of length 1 (o), 3 (□), 4 (A) 
and 5 (V). Simulation results are based on 1000 realizations 
of a chain of 20000 sites with J2/J1 > 1. 

Monte Carlo simulations of the T = Ising-Glauber 
chain yield results which are qualitatively consistent with 
these predictions. As shown in Fig. 1, the densities of the 
various elementary excitations decay at different tempo- 
ral rates. Over the last 2 decades of data, linear least- 
squares fits give 0.56 and 0.40 for the exponent associated 
with the density of l's and 3's, respectively. However, 
there is considerable curvature in the data and the linear 
fit is not indicative of the asymptotic behavior. In fact, 
the data for the density of l's is curved upward while 
that for the density of 3's is curved downward; these are 
suggestive of an asymptotic exponent of 1/2. However, 
a serious quantitative test of Eqs. ( |i"o| ) would require ex- 
tensive simulation, since the corrections to the local ex- 
ponents will vanish only as 1/ Int. We thus interpret the 
data as being consistent with (10). Similar least-squares 
fits give exponents for the density of 4's and 5's as 0.84 
and 1.42. In this case, both sets of data exhibit down- 
ward curvature. These two features are qualitatively con- 
sistent with the hypothesis that 4's and 5's are equivalent 
to [3, 3] and [3, 3, 3] resonances, respectively, so that their 
densities should scale as A 2 and A 3 /ln(l/A). 

To provide an additional insight into the basic nature 
of this coarsening process it is worthwhile to consider the 
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case J i = 0, where the chain breaks up into two inde- 
pendent antiferromagnetic sublattices. For each sublat- 
tice, the dynamics must coincide with the usual i -1 / 2 
Glauber coarsening. It is instructive to see how this be- 
havior arises within our picture of elementary excitations. 
The crucial feature for the case J\ = is that the process 
1 + 1 — * 3 + 3 now involves no energy change (see Eq. (||)). 
Therefore, the reverse process 3 + 3 — > 1 + 1 occurs freely. 
The presence of this additional process now gives the rate 
equations 



A = -3A 3 - 2A 2 
B = A 3 + 2A 2 



2B Z 
2B 2 



AB, 
AB. 



(11) 
(12) 



The cubic term turns out to be negligible and the result- 
ing asymptotic behavior is A(t) = B(t) ~ 1/t. Similarly, 
for the one-dimensional system, the density of l's and 3's 
become identical and both decay as t~ x l 2 . The presence 
of this additional reaction process and the corresponding 
change in the dominant terms in the rate equations shows 
clearly that the dynamics of the Ising-Glauber chain with 
competing interactions is in a different universality class 
than that of the purely antiferromagnetic spin chain. 



B. Weak Antiferromagnetism: J2 < Ji 

To complete the discussion of the Ising-Glauber chain 
with competing interactions, let us extend the consider- 
ations discussed above to J 2 < J\. For J 2 < Ji, the 
basic new dynamical feature is that nucleation of iso- 
lated domains in energetically forbidden at T — 0. Fur- 
ther, any l's in the initial population of elementary ex- 
citations quickly disappear through the annihilation pro- 
cesses 1 + 1^3 + 3 and 1 + 3^0. However 3's saturate 
at a finite density, as there is no mechanism for their re- 
moval. Resonances such as [3, 3], [3, 3, 3], etc., form freely 
by collisions between 3's. However, high-order resonances 
cannot decay by nucleation but merely re-fragment into 
lower-order resonances. 

In the case of J 2 < Ji/2, domain walls become the 
basic excitations, as in the ferromagnetic Ising-Glauber 
chain. However, domain walls repel since there is an 
energy cost AE = 4J 2 to bring two domain walls closer 
than two lattice spacings apart. Interestingly, for all J 2 < 
Ji, the equilibrium ground state cannot be reached via 
single spin-flip Glauber dynamics. For Ji/2 < J 2 < J\, 

the alternating . . . + H h + • • ■ ground state cannot be 

attained because there is no mechanism for domains of 
length > 3 to break up; for J 2 < Ji/2, the ferromagnetic 
ground state cannot be attained because of the repulsion 
of domain walls. 

The situations where J 2 = J\/2 and J 2 = J\ exhibit 
additional new dynamical features and each requires a 
separate treatment. In the former case, an arbitrary ar- 
ray of domains which all have lengths > 2 forms a ground 
state M. The Glauber single spin-flip dynamics now has 



the following geometrical manifestations. First, it is en- 
ergetically favorable for isolated 1-domains to disappear 
by "anti-nucleation" . Second, domain walls can diffuse 
freely as long as they are more than two lattice spacings 
apart. They repel each other with infinite strength at 
zero temperature, so they cannot annihilate. Most inter- 
esting is the dynamical behavior of strings of consecutive 
antiferromagnetic spins. In this case, the flipping of a sin- 
gle spin gives rise to a 3-domain. This can be viewed as 
the replacement of the triplet of antiferromagnetic spins 
by a ferromagnetic "trimer" . This replacement process 
of 3 spins by a trimer can continue inside an initially 
antiferromagnetic region until it is converted into a fer- 
romagnetic domain. 

This dynamics is therefore equivalent to the random 
sequential adsorption of trimers which can overlap. The 
dynamics of this adsorption process can be solved accord- 
ing to well-known procedures |l2j . Let P n (t) denote the 
probability that exactly n consecutive sites are occupied 
by 1-domains, that is, an antifcrromagnetically ordered 
string of length n. The probabilities P n (t) evolve accord- 
ing to the rate equations 



P n = -(n + 2)P n + 2 J2 P r 



(13) 



j=n+l 



The loss term arises because the adsorption of a trimer 
whose center coincides with any of the n sites of the string 
or the two sites adjacent to the string destroys the n- 
string. Similarly, the gain term arises from processes in 
which the adsorption of a trimer onto a larger string leads 
to the creation of an n-string. 

To solve the rate equations, note that the ansatz 



P n (t) = F(t)f n+2 (t) 



(14) 



transforms the infinite set of differential equations (|l3| ) 
into the pair of equations 



/ = -/, F 



2Fp_ 
1-/' 



(15) 



Consider now an initially disordered Ising chain (corre- 
sponding to a quench from T = 00 to T = 0). An anti- 
ferromagnetic string of length n occurs with probability 
P n (0) = 2~ n - 3 , which implies 



F(0) = /(0) = 



1 



Solving Eqs. ( |l5| ) subject to ( |l6| ) yields 
F(t)=2[ l-~e-" ] .-m> 



4 4 



(16) 

(17) 
(18) 



Therefore domain walls inside antiferromagnetic strings 
quickly disappear and the total density of domain walls 
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reaches a saturation level. Since the subsequent dynam- 
ics does not allow for processes which change the number 
of domain walls, the system continues to explore a sec- 
tor of the phase space that contains all ground states 
with the same number of domain walls. The system falls 
into its asymptotic sector exponentially in time. Different 
sectors are mutually disconnected, and the sector eventu- 
ally reached by the system depends on initial conditions. 
The dependence on initial conditions and exponential re- 
laxation are outcomes of the lack of ergodicity in the 
Glauber dynamics of the Ising model with competing in- 
teractions when J2 = Ji/2. 

Consider finally the marginal case J2 = J\. The cru- 
cial difference with the case of strong antiferromagnetism, 
•h > J%, is that the process 3 + 3 + 3 — > 1 now involves 
no energy gain and thus the reverse process 1 — > 3 + 3 + 3 
can also occur freely. The rate equations (||) and (^|) are 
thus modified to 

A~W-3R-r, (19) 
B ~ R-r-B. (20) 

To determine the asymptotic behavior, it proves conve- 
nient to transform these equations to 

A - B ~ -4(R - B), (21) 
A + 3B ~ -4r. (22) 

We now ignore B on the left-hand sides. Combining the 

estimate r ~ BA 2 with Eq. @ gives B A/ A 2 . This 

implies B —A, so that Eq. ( pl| ) now gives B ~ R, 
or A ~ — yl 5 /ln(l/A). We therefore conclude that when 
J 2 = J\ the density of elementary excitations are 

/In A 1/4 / 1 \ 1/4 

A(t)~( T ) and B( t )~( TO ) . (23) 

Thus in the marginal case of J2 — Ji, the Ising-Glauber 
chain still coarsens, but at a much slower rate than in 
the case of strong antiferromagnetism, J 2 > J\. 

III. THIRD- AND MORE DISTANT-NEIGHBOR 
INTERACTION 

Our general approach can naturally be applied to 
longer range antiferromagnetic interactions. We first out- 
line basic features of coarsening for an antiferromagnetic 
third-neighbor interaction; the behavior for arbitrary 
range antiferromagnetic interaction follows inductively. 
The Hamiltonian now is H = — J2i(Ji s i s i+i ~ ^3 s iSj+3) 
and for sufficiently strong antiferromagnetic interaction 
the ground state consists of alternating domains of length 
3, ... + + H h + + Starting from the ferromag- 
netic up state, the evolution to the ground state again 
proceeds by a two-stage process when J 3 > J\. While 
the ground state occurs when J3 > Ji/3, we shall em- 
ploy the stronger inequality J 3 > J\ in the following to 



guarantee that this ground state is accessible via single 
spin- flip Glauber dynamics. 

In the initial nucleation stage, it is energetically favor- 
able for a spin within an large up domain to flip if this 
spin is three or more lattice spacings from any domain 
boundary. Thus domains of length > 7 are unstable to 
such nucleation events. It is also energetically favorable 
for this isolated down spin domain to grow to size 3, as 
long as the expanding domain wall remains at least 3 lat- 
tice spacings from adjacent domain walls. At the end of 
the nucleation stage, therefore, all domains have length 
< 6. Further, domains of length 1 or 2 must be sur- 
rounded by domains of length 3; otherwise the central 
domain would expand until its size reached 3. 

A basic observation is that the true elementary excita- 
tions are domains of length 2 and of length 4, as these are 
the only objects which diffuse freely within a stable sea 
of 3's. All other defects are resonant states of these two 
elementary excitations. To determine the nature of the 
coarsening, first consider the resonances of 4-domains. 
For example, a 5-domain is formed when two 4's meet 
and interact, so that one 4 shrinks to length 3, while the 
other grows to length 5. There is no energy cost associ- 
ated with this process, so that a 5-domain can be viewed 
as a [4,4] resonance. Similarly, a 6-domain is a [5,4], or 
equivalently a [4, 4, 4] resonance. Finally, a 7-domain may 
be produced by the conversion of 6 + 4 — > 7 + 3; thus the 
7 is a [4, 4, 4, 4] resonance. At the center of the 7-domain, 
a single spin can flip, thereby nucleating a 1-domain and 
two surrounding stable 3 domains. The penultimate sto- 
ichiometry of this process is therefore, 4 + 4 + 4 + 4^ 1. 

Conversely, consider the resonances and interactions 
associated with 2-domains. When two 2's, within a sta- 
ble sea of 3's, meet, a 1-domain is formed as indicated 
by 

+ + H + H h + + 

+ + H 1 h + +. 

Since there is no energy cost associated with this pro- 
cess, a 1-domain is simply a [2, 2] resonance. When three 
2's meet, there are several possible zero-energy-cost out- 
comes. If the interior spin of the outer 2-domain flips, 
then the result is 2 + 2 + 2 1 + 3 + 2. Since the 1-domain 
is a [2, 2] resonance, this process can be considered as the 
first step in separating the three initial domains. How- 
ever, if one of the spins in the middle domain flips, then 
the outcome is 2 + 2 + 2^3 + 1 + 2. Once a 3, 1, 2 state 
is reached, it is energetically favorable for the central iso- 
lated spin to flip thus giving 3 + 1 + 2 — > 6, i.e., a [4, 4, 4] 
resonance. This last step is accompanied by the energy 
loss — 4Ji. The outcome of more than three 2's meet- 
ing can be obtained by grouping the 2's into triplets and 
analyzing the outcome of each triplet in series. Finally, 
when a 2 and a 4 meet, it is energetically favorable for 
2^3 and 4^3. This can be viewed as the two-species 
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annihilation 2 + 4 — > 0, since the two 3's formed in the 
reaction belong to the stable ground state. 

From these basic processes, the governing reactions for 
this system are 



4 + 4 + 4 + 4 -> 1 
2 + 4 -> 

2 + 2 + 2 -> 6 



2 + 2 



(24) 



For these reactions, the associated rate equations for the 
density of 4's (A) and 2's (B) are 



.4 
B 



AA 4 
2A 4 



3B 3 
3B 3 



AB, 
AB. 



(25) 
(26) 



The structure of these equations is similar to the second- 
neighbor interaction case. Following the same reasoning 
as used previously, we find A(t) ~ i -1 / 3 and B(t) ~ t~ l . 
Similarly, we may adapt the rate equation above to de- 
scribe the system in one dimension by following the ap- 
proach used to write Eqs. (ph and (S). This leads to 



.4 
B 



AA 4 



A 2 B. 



2A 4 - A B, 



(27) 
(28) 



with the asymptotic behavior A(t) ~ t 1 / 3 and B(t) 
t~ 2 ' 3 . 



of asymptotic behavior. We merely point out that our 
picture of elementary excitations allows one to express 
all densities in terms of the density of 4's. This predicts 
that the density of l's scales as A 4 , the densities of 2's 
and 5's scale as A 2 , and the density of 6's scales as A 3 . 
This is only marginally consistent with the data, a fea- 
ture which we attribute to slow approach to asymptotic 
behavior. 

Finally, we may apply a similar geometrical picture 
of elementary excitations to treat the general case of 
the Ising-Glauber chain with competing first- and n th - 
neighbor interactions when J„ > J\. In this case, the 
basic excitations are domains of n + 1 spins (A) and do- 
mains of n — 1 spins (B), each of which diffuses freely 
within a stable ground state sea of alternating ferromag- 
netic strings of length n. Other types of excitations are 
resonances of these elementary excitations. The basic ki- 
netic mechanisms that govern the coarsening of the spin 
system are the (n + l)-particle annihilation of contigu- 
ous groups of A excitations via (n + 1)A — > 1, and the 
two-species annihilation A + B — > 0. An analysis of the 
corresponding rate equations again indicates that these 
two processes are of the same order of magnitude. Con- 
sequently, the ground state is approached as t~ x l n and 
the densities of the basic excitations decay according to 
A{t) ~ t- x l n and B(t) ~ H 1 *- 1 )/". 
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FIG. 2. Density of domains of length 1 (*), 2 (o), 4 (□), 5 
(A), and 6 (V). Simulation results are based on 500 realiza- 
tions of a chain of 20000 sites with J3/J1 > 1. 

As in the case of the Ising chain with first- and second- 
neighbor interactions, simulations for the densities of the 
various elementary excitations indicate that they decay 
at different temporal rates (Fig. 2). Over the last 2 
decades of data, linear least-squares fits give 0.80 and 
0.25 for the exponent associated with the density of 2's 
and 4's, respectively. Once again, the sense of the cur- 
vature in these two data sets is consistent with the re- 
spective asymptotic exponents of 1/3 and 2/3. However, 
even more so than in Fig. 1 the linear fit is not indicative 
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We have investigated the coarsening kinetics of the 
Ising chain with single spin-flip dynamics when a more 
distant- neighbor antiferromagnetic interaction — J n com- 
petes with the nearest-neighbor ferromagnetic interac- 
tion Ji. For J n > Ji, this competition leads to slower 
zero-temperature coarsening compared to the case of 
nearest-neighbor ferromagnetic interactions only. Other 
types of non-universal relaxation phenomena have been 
reported for alternating interactions and other modifica- 
tions of the pure Ising chain |ll|. The case of the com- 
peting interaction, however, is amenable to an intuitively 
appealing description in terms of elementary excitations 
which makes clear the mechanism for the new types of 
coarsening kinetics. It is intriguing that the nature of the 
elementary excitations and the spectrum of resonances of 
the system are not obviously connected with the micro- 
scopic interaction of the spin system. 

For general n th -neighbor antiferromagnetic interac- 
tion, the elementary excitations are ferromagnetic strings 
of n+ 1 and n— 1 spins. The former interact and disappear 
through (n + l)-particle single-species annihilation while 
both excitations mutually annihilate when they meet. 
These two processes are of the same order of magnitude 
so that the rate of the overall coarsening process can be 
viewed as being limited by (n + l)-particle annihilation. 
This leads to a coarsening which proceeds as t~ x / n . The 
marginal case J n — J\ admits additional microscopic pro- 
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cesses which leads to even slower coarsening. Finally, it 
is interesting to note that Glauber dynamics has only 
tenuous connection with equilibrium properties of the 
system. That is, the dynamical criticality at J n = J\ 
is disconnected from the corresponding equilibrium be- 
havior, where fcrromagnctism occurs for J n < J\/n, a 
ground state of alternating domains of n ferromagnetic 
spins occurs for J„ > Ji/n, and an infinitely degenerate 
ground state consisting of alternating domains of > n 
spins occurs for J n = J\/n. 
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